Far From Threshold Buckling Analysis of Thin Films 
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Thin films buckle easily and form wrinkled states in regions of well defined size. The extent of 
a wrinkled region is typically assumed to reflect the zone of in-plane compressive stresses prior to 
buckling, but recent experiments on ultrathin sheets have shown that wrinkling patterns are signif- 
icantly longer and follow different scaling laws than those predicted by standard buckling theory. 
Here we focus on a simple setup to show the striking differences between near-threshold buckling 
and the analysis of wrinkle patterns in very thin films, which are typically far from threshold. 



The growing interest in developing technologies at 
smaller and smaller scales has posed new questions and 
challenges for scientists to understand the mechanical be- 
havior of tiny structures. Engineered films with thickness 
ranging from nano to microscales and designed for dif- 
ferent applications are among the ubiquitous examples 
of flexible structures that buckle under very small loads. 
More interestingly, these buckling instabilities usually de- 
velop into wrinkled patterns that form a dramatic display 
of the applied stress field [1, 2]. Wrinkles align perpendic- 
ularly to the compression direction, depicting the princi- 
pal lines of stress and providing through their geometry 
new tools for mechanical characterization. 

Buckling theory is regularly used to understand these 
patterns in macroscopic plates when the deformations 
are small perturbations of the initial flat state. How- 
ever, it has been known since Wagner [3, 4] that plates 
buckled under loads well in excess of those necessary to 
initiate buckling show an asymptotic state different from 
the one observed near threshold. The stress nearly van- 
ishes in the compression direction and the plate acquires 
fine wrinkles that mark the region where the compressive 
stress has collapsed. This asymptotic state is unusual in 
macroscopic plates, but very likely to happen in very thin 
films, since their threshold load values are very small. 

A better insight into this "collapsed" wrinkled state 
was provided by the recent discovery of scaling relations 
between wrinkle wavelength, film thickness and applied 
tension in stretched films of rectangular shape [5]. This 
theory, and later applications [6] , assumed that the wrin- 
kle length is determined by the in-plane compressive re- 
gion prior to buckling. However, recent experiments and 
theoretical work show that the length of wrinkles in very 
thin films is significantly larger than predictions based on 
the stress field near the onset of buckling [7, 8], and thus 
indicate that our conceptual understanding of the far- 
from-threshold wrinkled state is still lacking. In this Let- 
ter we present a novel analysis of the far-from-threshold 
limit and predict a new scaling law for the extent of the 



wrinkled region in very thin sheets. Details of this asymp- 
totic theory will be published elsewhere [9] . 




FIG. 1. The classical Lame configuration where a mismatch 
between the inner and outer stresses produces a compressive 
region for r < L. 



In order to study the essential differences between 
the near-threshold (NT) and far-from-threshold (FFT) 
regimes, we focus here on the circular configuration 
shown in Fig. 1 where an annular film of inner radius Ri n 
and outer radius R ou t is stretched differentially by radial 
forces per length T out at r = R out and T m > T out at 
r = R in . Similar geometries have been used to study the 
wrinkling pattern under different types of central loads, 
such as the impact of fast projectiles [10], the deadhe- 
sion and wrinkling of a thin sheet loaded at a point 
[11], and the wrinkling and folding of floating membranes 
[7, 12, 13]. However, Fig. 1 exhibits the simplest load dis- 
tribution leading to wrinkling with variable length and 
wavelength. It is a classical problem of linear elasticity 
[14] to obtain the radial and hoop stresses of the planar, 
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axially-symmetric state from radial force balance: 
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The radial and hoop strains are, respectively, e rr — 
^§p-, £ee — ^r, and u r (r) is the only nonvanishing dis- 
placement. Here, Y = Et is the stretching modulus, E 
is the Young modulus, t the film thickness, and v is 
the Poisson ratio. Focussing for simplicity on the limit 
Rout ~> Rin, the solution is: 

orr = (Tom + AT -pr ) ; &ee = (T out - AT ) (3) 

where AT = Tj„ — T ou t. This solution is usually named 
after Timoshenko, although it was first studied by Lame 
[14]. In this solution the radial stress is tensile every- 
where, but if the stress ratio T in /T out > 2 the hoop stress 
becomes compressive for R in <r <L NT . Here: 

L NT Rin \/Ti n /T ou t — 1 Ri Rin (Ti n /T ou t) ^ ■ (4) 

The existence of compression leads to the formation of a 
buckled state if the film is sufficiently thin that it becomes 
energetically favorable to relieve the compressive stress 
by bending. The Lame problem is thus characterized by 
two independent dimensionless groups: 
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where B = Et 2 / 12(1 — v 2 ) is the bending modulus. Buck- 
ling is expected for a fixed value of r > 2 when e is smaller 
than a critical value e c . Our analysis assumes reducing 
e from NT (e < e c ) to FFT conditions (e < 1) by, for 
example, varying the film thickness t. 

Standard buckling theory in the NT regime consists of 
a stability analysis of the I s * Foppl-Von Karman (FvK) 
equation [4]: 
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where the out-of-plane displacement is assumed to be of 
the form ( = f(r) cos(m#), and a rr ,ogg assume their 
Lame form (3). While the energy and wrinkled extent in 
the NT regime are determined by Eqs. (3,4), bifurcation 
analysis (assuming f(r) is infinitesimal) yields the nu- 
merical value of e c and the exact number of wrinkles to 
of the emerging pattern [15]. Simple dimensional anal- 
ysis implies the general result to = g(r,e), where g is 
some dimensionless function. Near threshold, the Lame 
solution assures that a rr ~ ogg . If T out <C Ti n we expect 
m » 1 and one notices that the out-of-plane stretching 
force in the hoop direction is much larger than in the 
radial direction. Thus, the bending forces on the LHS of 



Eq. (6) are balanced by the out-of-plane hoop stretching 
on the RHS, and one obtains the NT scaling law: 



to - g(r)Rin yjT out /B , 
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for some g(r). This linear relation between the number 
of wrinkles and hole radius was confirmed in Ref. [16] in 
experiments with macroscopic rubber membranes, which 
are expected to be under NT conditions. Similarly, 
Refs. [6, 16] reported that for large samples the wrinkle 
extent L increases as the square root of the inner tension, 
in agreement with Eq. (4). 

Let us turn now to our main focus: the FFT regime e<§; 
1. Motivated by experiments [7] and following Wagner's 
ideas, we assume that the sheet is composed of two parts: 
a wrinkled region in Ri n < r <L with a collapsed hoop 
stress ogg — > 0, and an outer annulus L < r < R out in 
which the sheet remains planar with stresses following the 
Lame form (3) appropriately modified. We shall prove 
below that for e -c 1 a state with L > R in is energetically 
favorable to the Lame state, which corresponds to L = 
R in - Thus, wrinkling is a mechanism for releasing elastic 
energy in the film. For R in <r<L, Eqs. (1,2) yield: 

a rr = T in —^- ; e rr = (T in /Y)—^- ; egg = —ve rr . (8) 
r r 

Since the radial stress must be continuous we find that 
the stresses in the outer annulus have the Lame form (3) 
with Rin^L and AT->AT(L) = T in R in /L — T out . For 
a given wrinkle extent L, the FFT stresses are now fully 
characterized by Eqs. (3,8). Moreover, the radial dis- 
placement u r (L) must also be continuous, and since the 
Lame solution at r > L implies a link between u r (L) and 
a rr (L) one obtains the radial displacement at r<L [9] : 

u r (r) = Rin (Ti n /Y) log + u r (L) , with 
u r (L) = [2L (T out /Y) - (1 + v)R in (T m /Y)}. (9) 

We shall determine the actual wrinkle extent by mini- 
mizing the energy over all allowed values of L. 

Before turning to energy calculations, let us highlight 
some important aspects of the FFT solution. First, as 
Eq. (8) indicates, there is a pure traction along the ra- 
dial direction producing a Poisson effect egg = —ve Tr < 
in the azimuthal direction and reducing the perimeter at 
radius r by a total length — 2-KrvTi n /Y . This local con- 
traction shows that for v ^ the film is not inextensible 
in the azimuthal direction as was assumed in [5]. How- 
ever, a similar constraint arises: there is an excess in 
length when this contraction is not compatible with the 
geometrical shortening of the perimeter length 2iru r (r) 
that is generated by the inwardly radial displacement 
u r < 0. In order to reduce stretching energy, this ex- 
cess of length is relieved by out-of-plane displacement 
that is highly oscillatory in the azimuthal direction. Us- 
ing the weakly nonlinear strain-displacement relations [4] 
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Assuming again that C = f( r ) cos(m(9) one finds: 
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The positivity of the LHS of Eq. (11) manifests that out 
of plane deformation is possible only if there is an excess 
of length of the circle perimeter. An algebraic manipu- 
lation of Eqs. (9,11) shows that the RHS is positive for 
r < L only if L<Ri n Ti n /2T out , thus providing an upper 
bound for the wrinkled extent. Additionally, Eq. (11) 
indicates that the product mf(r) must remain finite in 
the FFT limit, in contrast to the NT regime where f(r) 
is infinitesimal. Indeed, our FFT analysis is based on 
an asymptotic series in which the expansion parameter 
is y/e, in contrast to the NT regime, where the small 
parameter is the distance to threshold e — e c [9]. 

In order to determine the wrinkled extent, L, we com- 
pute the elastic energy Ue of the FFT stress field, where 

Ue = \ J dA (°Vr£rr + oeesee) ■ (12) 

A straightforward calculation using Eqs. (3,8) yields: 

U E = (tt/Y) {(1 - v){R out T out ) 2 + [R m T in f \n(L/R m ) 
+ 2(LT out - R ln T ln f - (1 - V )(R m T ln f + 0{l/R 2 out )} 

(13) 

Note that this energy does not include the costs of bend- 
ing and out-of-plane stretching of the membrane. These 
may be shown to be higher order contributions in e in 
the FFT limit [9], but are nevertheless crucial for ob- 
taining the exact limit of the number of wrinkles, m. 
Since our problem involves applied constant forces at 
r = Ri n , R ou t , we must minimize the mechanical en- 
ergy U = U E -W, where W = 2-K[T out R out u r (R out ) ~ 
Ti n Rin u r (Rm)] is the exerted work. Minimizing U as 
a function of wrinkle extent is analogous to fracture me- 
chanics problems, in which one minimizes the mechanical 
energy as a function of crack length under constant load 
conditions [17]. Like cracks, wrinkles provide a route 
for the release of elastic energy. Using a general result 
from elasticity theory for bodies under constant external 
loads [17], we find that W — 2Ue, and hence U = —Ue- 
In order to minimize U we notice that its first derivative 
d L U = -Tt(2LT out - R m T m ) 2 /(YL) < is zero for: 



-^FFT 



2 -7oi. 



(14) 



Interestingly, this result is identical to the upper bound 
that we found above, and also assures continuity of the 



hoop stress aeg(L FFT ). Thus, energy minimization natu- 
rally yields a value for the stress at the tip of the wrin- 
kles that smoothly matches the flat region in the film to 
the highly wrinkled one. The quadratic form of the first 
derivative and Fig. 2 show that Eq. (14) corresponds to 
an inflection point of U , but the upper bound guarantees 
that this is the actual length in the FFT limit. Fig. 2 
also indicates that at the FFT limit, where corrections 
to the mechanical energy U are ~ -^/e, the energy at the 
inflection point U(L FFT ) is lower than the Lame value 
U(Rin)- Thus the energy is lowered from the NT (Lame) 
value to the FFT asymptotic limit, as e is reduced from 
its critical value e c to e<Cl. Equation (14) is the central 
result of our analysis. It reflects a linear scaling of the 
extent of the wrinkles in the FFT regime with the ra- 
tio Ti n /T ou t, in sharp contrast to the square root scaling 
that characterizes the NT limit, Eq. (4). 
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FIG. 2. The mechanical energy U in the FFT limit (nor- 
malized by Y / (RinTout) 2 ) as a function of the wrinkle ex- 
tent L (normalized by Ri„). Here we took T in /T ou t = 8, 
Rout = lORin and v — 1/3. An inflection point exists at 
L FFT = Ti n /2T ou t — 4; the upper bound implies that this is 
the actual extent of the wrinkles. When reducing e from e c 
in the NT to the FFT regime (e^C 1) the energy decreases by 
an amount AU from the upper dashed line (the energy of the 
Lame solution) to the value at the inflection point, [/(Lfft). 

Let us turn now to the number of wrinkles m in the 
FFT limit e ->• 0. While Eq. (11) reveals that the prod- 
uct m/(r) must remain finite, we do expect |/| to become 
smaller and hence m diverges as the FFT regime e <C 1 is 
approached. This agrees with experimental observations 
[5, 7], and is also a direct consequence of our theory [9]. 
Physically, this is correlated with the fact that the dom- 
inant out-of-plane forces in the FFT regime are bending 
and compression in the azimuthal direction and stretch- 
ing in the radial direction. To leading order in e, the I s * 
FvK Eq. (6) in the FFT regime is [9]: 
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Balance of these three forces implies the hoop stress scal- 
~ V BT out I Rin in contrast to the NT regime, 



mg a e , 
where 



is a function of the exerted stresses T™,T, 
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but is independent B, see Eq. (3). One should notice that 
this result can be expressed as <Jee/T out <~ e 1 ^ 2 , and re- 
flects the collapse of compressive stress in the FFT regime 
that is the basis of our analysis. The same balance im- 
plies the FFT scaling law for the wrinkles number 

m ~ k(r)(Rl T^/B) 1 ' 4 ~ e" 1 / 4 (16) 

for some fc(r), in sharp contrast to the NT scaling (7). 
This scaling law has already been predicted in [5] , and is 
strongly supported by experimental observations on ultra 
thin-sheets, where e is estimated to be below 10 -6 [7]. 

In the experiments of [7], a very thin circular sheet is 
floating on water, subject to surface tension T out = 7 
at its perimeter r = R ou t- A liquid drop is placed at 
the center, deforming the sheet beneath it, and exerting 
an in-plane tensile force T„ at the contact line r = Ri n . 
The problem is thus analogous to the Lame problem, but 
the determination of T in is a subtle problem that has not 
been resolved yet [8]. Although our result (14) cannot 
be directly compared to the experiments of [7], it does 
provide a new answer to a puzzle raised by the empir- 
ical rule found there: L FFT w Cl\/Y j^R in with Cl a 
numerical constant. The authors of [7] assumed the NT 
scaling, Eq. (4), and concluded that their results indicate 
that Ti n is "independent of surface tension, which is im- 
plausible". However, the FFT scaling, Eq. (14), shows 
that the empirical law is consistent with Ti n <~ y/Yj, 
suggesting instead that the in-plane tension exerted by 
the drop at the contact line is affected "equally" by the 
surface tension and the stretching modulus. 

We conclude by making several comments on the na- 
ture of the FFT analysis that we have developed here. 
First, similarly to the NT analysis, we note that the wrin- 
kle extent (14) is determined by minimizing the in-plane 
stretching energy which is much larger that the bend- 
ing and out-of-plane stretching energies that determine 
the number of wrinkles (16). This leads to the predic- 
tion that in wrinkling patterns the length is more robust 
than the number of wrinkles. Second, our theory of the 
FFT regime (e -C 1) is similar to a "membrane limit" 
analysis which considers an elastic sheet with vanishing 
bending modulus (i.e. e = 0) [4]. Our analysis, however, 
clarifies the singular nature of this limit: while the NT 
state can be obtained by regular expansion around the 
planar Lame state [15], the FFT state is obtained by an 
asymptotic expansion around a limit in which to diverges. 
Finally, one may wonder whether our theory enables an 
exact computation of the asymptotic number of wrinkles 
to beyond the scaling law (16). In order to do this, one 
must use Eqs. (8,9,11) to compute the bending energy 
Ub = (B/4) drm 4 f 2 (r) and the out-of-plane stretch- 
ing energy Us = (1/4) J^, dr a rr (r)f \r) 2 as a function 
of to, and determine m as the minimizer of these sub- 
dominant energies. One finds, however, a divergence of 
f'(r) and Us as L — > L FFT for all m. Initially, this obser- 



vation may lead one to doubt the validity of our results. 
However, careful thinking reveals the source of this di- 
vergence [9]: our matching conditions at r = L assume 
a direct transition from the planar state at r > L to the 
fully-collapsed region at r < L. The divergence of Us is 
cured if one replaces the "pointwise" matching by a nar- 
row boundary layer near r w L in which the two regimes 
are smoothly matched. The energetic cost of this bound- 
ary layer is essential to find to, but requires a nontrivial 
calculation which is currently underway. 

To summarize, we introduced an analysis of thin sheets 
far from buckling threshold. The central assumption of 
our theory, which should also be valid for more compli- 
cated set-ups, is the collapse of compressive stress in the 
wrinkled region. For the simple Lame geometry we found 
a new scaling law for the wrinkle length and showed that 
it explains a puzzle arising from previous experiments [7] . 
We anticipate that our work will lead to understanding 
of far-from-threshold wrinkling patterns under many dif- 
ferent load configurations. 

We acknowledge support by the Petroleum Research 
Fund of ACS (B.D.) and NSF-MRSEC on Polymers at 
UMass (R.S.). D.V. is supported by an Oppenheimer 
Early Career Fellowship. M.A.B. and E.C. acknowledge 
the support of CNRS-Conicyt 2008. E.C. thanks Fondc- 
cyt project 1095112 and Anillo Act 95. 



[1] N. Bowden, S. Brittain, A. G. Evans, J. W. Hutchinson, 
and G. M. Whitesides, Nature 385, 149 (1997). 

[2] J. Genzer and J. Groenewold, Soft Matter 2, 310 (2006). 

[3] H. Wagner, Z. Flugtechn. Motorluftschiffahrt 20, Nos. 
8-12 (1929). 

[4] E.H. Mansfield, The Bending and Stretching of Plates, 

(Cambridge University Press, Cambridge, 1989). 
[5] E. Cerda and L. Mahadevan, Phys. Rev. Lett. 90, 074302 
(2003); E. Cerda, K. Ravi-Chandar, and L. Mahadevan, 
Nature 419, 579 (2002). 
[6] E. Cerda, J. Biomechanics 38, 1598 (2005). 
[7] J. Huang et al. Science 317, 650 (2007). 
[8] D. Vella, M. Adda-Bedia and E. Cerda, Soft Matter 

(2010, accepted). 
[9] B. Davidovitch and E. Cerda (in preparation). 
[10] R. Vermorel, N. Vandenberghe, and E. Villermaux, Proc. 

R. Soc. Lond. A 465, 823 (2009). 
[11] J. Chopin, D. Vella, and A. Boudaoud, Proc. R. Soc. 

Lond. A 464, 2887 (2008). 
[12] J. Huang, N. Menon, and T. P. Russell, Soft Matter 

(2010), submitted. 
[13] D. P. Holmes and A. J. Crosby, Phys. Rev. Lett. 105, 
038303 (2010). 

[14] S. P. Timoshenko and J. N. Goodier, Theory of Elasticity, 
McGraw Hill, 1970. 

[15] G. G. Adams, J. Applied Mech. 60, 520 (1993). 

[16] J.-C. Geminard, R. Bernal, and F. Melo, Eur. Phys. J. 
E 15, 117 (2004). 

[17] B. Lawn, Fracture of Brittle Solids (Cambridge Univer- 
sity Press, 2nd edn., Cambridge, 2004). 



